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DISCLAEMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
COMMON, as to the accuracy and functioning of the program 
and related program material, nor shall the fact of distribution 
constitute any such warranty, and no responsibility is assumed 
by the contributor or COMMON, in connection therewith. 
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COMMON USER3 GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter, ink or pencil) 

^frogram No. Date 

Program Name: 

1. Does the abstract adequately describe what the program is and what Yes No_ 

it does ? 

Comment ^ 

2. Does the program do what the abstract says 7^ Yes No_ 

Comment 

3. Is the description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No_ 

C om m ent 

Are the Sense Switch options adequately described (if applicable)? Yes No__ 

Are the mnemonic labels identified or sufficiently understandable? Yes No_ 

Comment 



CI 



5. Does the source program compile satisfactorily (i£ applicable)? Yes No_ 

Comment 



; Does the object program run satisfactorily? Yes No_ 

Comment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Comment \ 

8. Does the Program meet the minimal standards of COMMON? Yes No_ 

Comment 



9. Were all necessary parts of the program received? Yes 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to; Your Name 

Company 



NO' 



Mr. Richard L. Pratt 

Data Corporation Address 
7500 Old Xenia Pike 



^ Dayton, Ohio 45432 Users Group Code 

THIS REVIEW FORM IS PART OF THE COMMON ORGANIZATION'S PROGRAM REVIEW ANT 
EVALUATION PROCEDURE. N0NMEMBER3 ARE CORDIALLY INVITED TO PARTICIPATE 
IN THIS EVALUATION. 
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DOES 

(Disk Oriented Equation Solver) 



Author: Mr. Robert Stephenson 

University Computing Center 
University of Cincinnati 

Date : 1/29/65 



Membership Code: 30l5 



"Modifications or revisions to this program, as they 
occur, will be announced in the appropriate Catalog of 
Programs for IBM Data Processing Systems. When such 
an annovmcement occurs, users should order a complete 
new program from the Program Information Department." 
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PROGRAM ABSTRACT 

TITLE (If subroutine, state in Title): "DOES" (Disk Oriented Equation Solver) 
Subject Class if ication: 5«o 



2. Author; Organization: Mr. Bnh»rt Rti>ph>wBnn 



Pniversity Compttting Center o f the Jniversitv cf Cincinnati 

Date: 1/28/65 ^Users Group Membership Code: 3015 

3. Direct Inquiries to Name: pr. Carl F. Evert. Director 

University of Cincinnati Phone: 75-2 333 

4. Description/Purpose: (5. Method; 6. Restriction/Range; When Applicable): 

This proKTam utilizes the I3II Disk to extend the capacity 

of the 1620 for the solution of simultaneous algebraic equations or 

matrix in Tersion. 200 x 200 matrices may be inverted. 



7. Specifications (Check or fill in appropriat'? spaces): 
a. Storage used by progrsim: ^ok and 1311 Disk 



b. Equipment required by program: Card X ; Magnetic Tape ; Niimber of 

Drives ; Paper Tape ; Disk Pile x ; Number of Drives 1 ; TNS, 

TNF, MF X ; Auto divide x ; Indirect addressing x ; Floating Point 

Hardware x ; 1620 Model I x ; Model n ; 1443 Printer ; Index 

Registers ; Binary Capabilitit'i; ; Other (specify) 



Can program be used on lesser macni:i' ; Yea • Specify which requirements 
ran be easily removed "Additional p*mark«." . 

c. Programmed in: Fortran without Forn.at ; Fortran with Format ; 

Fortran II ; Ot.her Fortran (specify) ; SPS (specify assembler used) 

sps-iTD : Other (specify) 

d. Type of Program.: MairJine, complet e X ; Subroutine ; If subroutine, for 

use with SPS (specify type of SPS) ; .'-"ortran (specify Vj^^e of Fortran 

; Other (specify) 

Additional Remarks: Program could be used on a 20K machiae with no difficulty 
and with reasaemble^ it could be used on machines without the automatic 

floating pniaf; f>a».»n»»- , 
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PROGRAM WRITE-UP 



a) DOES (Disk Oriented Equation Solver) 

b) Date: 1/29/65 

c) Dr. Carl P. Evert, Director; University of Cincinnati; 
1+75-2333; Membership Code - 3015 

d) Program Description 

The technique used for this inversion is based on a process 
described by Sherman.^ 

Sherman states : Given a matrix [k] differing from another 
matrix (jal only in the elements of one column, the inverse of [a3 , 
P3» can be found knowing the inverse of (jaj , [b] , and the elements 
of the kth column of TaJ . Where the kth column is the column of 
difference betweenf AJ and [&"] . The expressions for calculating 
the elements of [BJ are: 



i = 1,2,3. 



= J. ^nr^rk 
r=l 



n = 1,2,3. 



,,N, 



Bij = bij - ZiBkj 



l,2..k-l,k+l,...N. 
1,2,3,. N. 



N = 



Bab= 

\b= 

^rk= 
Z_ = 



The size of the matrices. 

The element of [b] located in the ath row and the bth column. 
The element of [bj located in the ath row and the bth column. 
All the elements in the kth column of [A^ , since r=l,2, 3, . . .N. 
The matrix product of the nth row of [b^ and the kth column of (^A]. 
Consider the following example: 



10 
10 
1 



k = 1 and N = 3. 



10 
10 
1 



[A] = 



o~l 

1 

ij 



Sherman, Jack; "Adjustment of an Inverse Matrix Corresponding to 
a Change in the Elements of a Given Column," U.S. Department of 
Commerce, National Bureau of Standards, Applied Mathematics 
Series # 29. 
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z, = 



Z-, = 



[l o] . 
[o i] . 



= 7. 



= 8. 



Zp =[0 1 o]. 



= 5. 



«11= 


^ll/^l = 1/7- 


^12= ^12/^1 


= 0/7 


^21= 


bgi- Z^B^^ = 


- 5(1/7) = 


-5/7. 


B22= 


^22" ^2^2 ^ 


- 5(0) = 


1. 


B23= 


b23'" ^2^13 ~ ^ 


- 5(0) = 


0. 


^31^ 


^31" ^3^11 " ° 


- 8(1/7) = 


-8/7. 


832= 


b32" '^^^2.2. ~ ^ 


- 8(0) = 


0. 


B33= 


^33" ^3^13 " 


- 8(0) = 


1. 



^13' "13^ 



, , ^ r 1/7 o") 



The general Inversion program starts by generating a N x N 
identity matrix and storing it on the disk. The matrix on the 
disk is the ^h] matrix. The first column of data is read into the 
memory. An inversion process is then performed letting [a] be a 
matrix different from an identity matrix only in the first column, 
and that column is the one just read. As the inversion is 
performed, the |^B"] matrix is stored on the disk replacing the 
identity matrix. A second column of data is read into memory and 
(a^ matrix is now different from an identity matrix in the first 
two columns, however, it differs from a matrix Cal » whose 
inverse is now on the disk, only in the second column. The [b3 
matrix from the previous inversion now becomes the CbT matrix for 
this inversion. This process continues iintil all the columns of 
data have been read, at that time the matrix on the disk is the 
desired inverse of the input matrix. 

The common method of attach on problems of systems of 
simultaneous linear equations is to use a matrix inversion technique, 
The system 

^-ll^l"^ ^12^1"*" "•' *ln^ ^ °1 
a23_x^+ a22X2+ .... + a.^^.^ = 0^ 



nn 1 



■n2^2'^ 



can be represented in matrix form by 





^1 




"°1 




^2 
















X 

n 




c 

n 



•"ll ''12 ^In 
^21 *22 ••• *2n 

a , a ^ ... a 

nl n2 nn 

If we call the coefficient matrix A, the variable or unknown 
matrix x and the constant matrix c, we have 

AX = C. 

If we premultiply both sides by A"-'-, the inverse of A, we have 

A-l(AX) = A-lC 
or 

(A"^A)X = A"^G 
X = A~^C 

The problem is then to invert the Z matrix and multiply it 
into the C matrix. 



e ) Input/Output Formats 

The first card contains 1 parameter, the number of equations, 
punched in the first four columns and right Justified. * 

The next set of cards is the A or coefficient Matrix punched 
by columns . The format is like the Fortran format 5{Elij.,7,2X) , 
but unlike Fortran no deviation is allowed. Each new row is to 
be started on a new card, and all numbers are to be normalized. 

If desired, the program will compute the equation solutions. 
The G matrix is put on cards in the same manner as a column of the 
A matrix, and it follows the A matrix in input. 





/ 

A-matrix 


Parameter 
Card 








* WARNING: If this parameter 
is greater than 200 or is 
not right justified in the 
colvunns indicated, it is 
possible for the Monitor 
System to be destroyed. 



Output Format 



The incomplete and complete inverse are punched in the same 
format, 5(Eli|..7,2X) , as the input. The complete inverse has 
headers between the rows, but the incomplete inverse does not. 

The solutions are punched one to a card with the number and 
solution of each unknown. 



f ) Restrictions 

This program is designed to handle up to 200 simultaneous 
equations using an IBM 1620 - I3II system. The program was 
assembled on a 1620 with l\.0&. storage, automatic floating point 
hardware, and indirect addressing using SPS-IID, It could be 
used on a 20K machine with no difficulty and with reassembly, it 
could be used on machines without the automatic floating point 
feature . 

This program is compatible with Monitor I as it is issued 
from IBM. The first 20 cylinders of disk storage are used for 
the 200 equation problem, but for smaller problems, the cylinder 
usage can be coBQ>uted by 

NiMiber of Sectors = (Number of Equations )^ 
and there are 200 sectors per cylinder. 



g) Operating Instructions 

1. Cold - Start Monitor I 

2. Load Program and input cards into read stacker. 

3. Ready Punch 

\\., Push reader start 



Details 

Switch Settings 

At least one progrsun switch must be on to get any output. 

SWITCH 1 - On to make intermediate exit 
SWITCH 2 - On to compute and punch equation solutions. 
SWITCH 3 - On to restart after intermediate exit. 
SWITCH ij. - On to punch completed inverse matrix. 

All Check Switches to STOP 
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Intermediate Exit 

If desired, the program can be stopped during a run and the 
partially inverted matrix punched for restart later. To do this 
turn on Switch 1 at the desired stopping point. 

To restart place the incomplete inverse deck in the reader 
followed by the part of the A matrix not already read by the pro 
gram. If in doubt, the first card of the incomplete inverse 
should be printed. It contains 3 parameters, the number of 
equations, the next column to be read, and the number of sectors 
per column (always 20), 
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Flow Chart 



Clear 10 Areas 
Set Group Marks 
At Ends of 
Arrays 



Read 

Parameter 
Card 



Switch 
3 



Off 



fOn 



Store JJJ,NN,N 
Read Partial 
Inverse and 
Store on Disk 



K=JJJ 



Compute Sector 
Address of K^^ 
Row and Seek 



Store N 
Set JJJ=1 
N=20 

Record Identity 
Matrix on D isk 
I 



- ^+12 



Read Kth 
Column into 
Array A 



Read k'^" Row 
of Inverse into 
Array C 



Compute 
n 

Z=I Aid 
i 




-7- 



J=l 



c+c/z 
J=J+1 



<0 



J-K 



>0 



Write C 


Array 


As New Rth 


Row 


OD Disk 


and 


Seek 


Row 1 














L= 


1 





To 
P26-24 





o 



a 



a 



1=1 

IE 



Head Sow from 
Disk into Array 
B ana Seek Kext 
Row 




Off 



Punch Row of 
Inverse 




Switch 

2 



On 



Off 



°I=i ^III^III 



1=1+1 



GZ) 




Off 



Punch Equation 
Solutions 



c 



Halt and Branch 
to Monitor 



Oinin»r)i<;K qrifntpo fiOUATION <;OLVER 



01020 START 


TFM 


*+18 ♦ lO-l ♦ » INI T I AL I ZAT ION 


02402 


16 


02420 


-5502 




01 030 


TDM 




02414 


1 5 


00000 


00000 




1 040 


AM 


*-6»l 


02426 


1 1 


02420 


-0001 




01050 


CM 


»-18» 10+158 


02438 


14 


02420 


-5661 




01060 


BNH 


*-36 


02^50 


47 


02414 


01 100 




01061 


TFM 


•+18.A-9 


02462 


16 


02480 


-5688 




01062 


TDM 




02474 


15 


00000 


00000 




01063 


AM 


»-6,l 


02486 


11 


02480 


-0001 




01064 


CM 


♦-18.C+1990 


02498 


14 


02480 


J1691 




01065 


BNH 


»-36 


02510 


47 


02474 


01100 




- Ol_06^ 


TD . 


_B-11»QM_ 


02522 


25 


07688 


05687 




01067 


TD 


C-11»GM 


02534 


25 


09690 


05687 




01068 


TD 


0NE-11»GM 


02546 


25 


11692 


05687 




01070 


RNCP 


I0-1»»»READ PARAMETER CARD 


02558 


36 


05502 


00500 




01080 


SF 


lO-l 


02570 


32 


05502 


00000 




01090 


BNC3 


PI 


02582 


47 


02806 


00300 




-OllOQ P450_ 


SF 


I0+4.»»C0NTINUATI0N OF OLD JOB 


02594 


32 


05507 


00000 




nilin 


«;f 


10+7 


02606 


32 


05510 


oonon 




01120 


TF 


JJJ»I0+2»» JJJ IS COLUMN TO BE 


READ NEXT 02618 


26 


05666 


05505 




01125 


TFM 


JJJ-3»0.10 


02630 


16 


05663 


000-0 





01130 
01140 



01150 
01160 

01170 D0475 



TFM NN»20»9. NN IS SECTOR COUNT 02642 16 05494 00-20 

TF N. 10+10 .» N IS NUMBER OF ROWS AND COLUMNS IN MATRIX 

. 02654 26 05671 05513 

TFM MATl.O 02666 16 05491 -0000 

SK DCF 02678 34 05486 00701 

TFM MMtln READ PARTIAL INVERSE AND 02690 16 05676 -0001 



01180 
01190 
01200 



01205 
01206 
01210 



BTM 
TFM 
_WDGN 
A 

SK 
AM 



READ.B. ♦ 
CORE.B-9 
DCF 

mati/nn 

DCF 
MM,1 



RECORD ON DISK 



01220 C MM»N 
01230 P475 RNH D0475+12 
0124 B7 P9 _ 

01250 PI TF N»I0+2» .START OF NEW JOB 

02020 TFM NN,20,9 

02030 TFM JJJ«1 



02702 17 11858 
02714 16 05499 
02726 38 05486 
02738 21 05491 
02750 34 05486 
02762 11 05676 



-7699 
-7690 
00700 
05494 
00701 
-0001 



02040 

02045 

02050__ 

02060 

02070 

02080 



02774 24 05676 
02786 47 02702 
02798 49 03150 
02806 26 05671 
02818 16 05494 
Q283Q 16 05666 



05671 
01 100 
00000 
05505 
00-20 
-OQQl 



TFM 
TFM 

SK _ 
TFL 
TFM 
MM 



MATltO 
C0RE»S-9 

DCF 

B»ONE»« 
L»2 

L .10 



RECORD IDENTITY MATRIX ON DISK, 



02090 



SF 



02842 16 05491 
02854 16 05499 
02866 34 05486 
02878 06 07699 
02890 16 05686 
02902 13 05686 



95 



02914 32 00095 



-0000 
-7690 
00701 
11703 
-0002 
-0010 
00000 
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? 1 C r: 




AM 


99,B-10 




02926 


1 1 


00099 


-7689 








Th"L 


-99.ZFR0 




02938 


06 


0009R 


11713 








AM 


L»l 




02950 


1 1 


05686 


-0001 




0213'^ 




C 


L ,N 




02962 


24 


05686 


05671 




0? 1 AO 




BNH 


♦-72 




02974 


47 


02902 


01 100 




02150 




TFM 


M» 1 




02986 


16 


05681 


-0001 




0216O 




87 


P6 




02998, 


49 


03078 


00000 




021 




MM 


M.IO 




03006 


13 


05681 


-0010 




021B'^ 




<^F 


95 




.. 3018 


3? 


.00095 


..Q.0000_. 


_ . 


02190 




AM 


99»S-10 




30 30 


11 


00099 


-7689 




022 00 




TFL 


-99 »ONE 




03042 


06 


0009R 


11703 




02210 




SM 


99. 10 




03054 


12 


00099 


-OOlO 




022?o 




TFL 


-99. ZERO 




03066 


06 


0009R 


11713 




02230 


P6 


WDGN 


DCF 




03078 


38 


05486 


00700 








A 


MATl .NN 




03090 


21 


05491 


05494 




02245 




5K 


DCF 




03102 


34 


05486 


00701 




02250 




AM 


M»,l 




03114 


1 1 


05681 


-0001 




02260 




C 


M.N 




031?6 


24 


05681 


05671 




02270 




BNH 


»-132 




03138 


47 


03006 


01100 




02275 


P9 


TF 


K.JJJ.. 


BEGINNING OF MAIN LOOP 


03150 


26 


11718 


06666 




02280 




TF 


FAC .K.. 


COMPUTE SECTOR ADDRESS OF 


03162 


26 


11728 


11718 




02290 




SM 


FAC.l.. 


K TH ROW OF INVERSE 


03174 


12 


11728 


-0001 




02300 




M 


FAC.NN 




03186 


23 


11728 


05494 




03C1C 




5F 


97 




03198 


32 


00097 


00000 




03030 




TF 


MATl .99 




03210 


26 


05491 


00099 




03C40 




SK 


DCF 




03222 


34 


05486 


00701 




03050 




BTM 


READ. A.. 


RFAD K TH COLUMN OF MATRIX 


03234 


17 


11859 


-5697 




03060 




TFM 


CORE.C-9 




03246 


16 


05499 


-9692 




03070 




RDGN 


DCF 




03258 


36 


05486 


00700 




03090 




TFL 


Z.ZERO 




03270 


06 


11738 


11713 




03100 


D013 


TFM 


MM. 1 . . 


COMPUTE Z 


03282 


16 


05676 


-0001 




03110 




MM 


MM.IO 




03294 


13 


05676 


-0010 




03120 




SF 


95 




03306 


3? 


00095 


00000 




031 30 




AM 


99.C-10 




03318 


11 


00099 


-9691 




03 1 40 




TFL 


FAC .-99 




03330 


06 


11728 


0009R 




03150 




AM 


99.A-C 




03342 


11 


00099 


-400M 




03170 




FMUL 


FAC .-99 




03354 


03 


11728 


0009R 




03180 




FflOD 


Z .FAC 




03366 


01 


11738 


1 1728 




03190 




AM 


MM.l 




03379 


U 


05676 


-0001 




03200 




C 


MM.K 




03390 


24 


05676 


11718 




03210 


P13 


BNH 


D013+12 




03402 


47 


03294 


01 100 




03220 


0014 


TFM 


J»l ». 


DIVIDE ELEMENTS OF K TH ROW 


03414 


16 


11748 


-0001 




03230 




MM 


J .1 0. . 


OF INVERSE BY Z 


03426 


13 


11748 


-0010 




03240 




SF 


95 




03438 


32 


00095 


00000 




03250 




AM 


99.C-10 




03450 


U 


00099 


-9^,91 




0326O 




FDIV 


-99. Z 


12 


03462 


09 


0009R 


11738 





03270 




AM 


J.l 




03474 


11 


11748 


-0001 


03280 




C 


J.K 




03486 


24 


11748 


11718 


03290 


P14 


BNH 


D014+12 




03498 


47 


03426 


01100 


03300 




TFM 


CORE.C-9 


..RECORD NEW K TH ROW ON DISK 


03510 


16 


05499 


-9692 


04010 




WDGN 


DCF 




03522 


38 


05486 


00700 


04020 




TFM 


MATl .0 




03534 


16 


05491 


-0000 


04030 




SK 


DCF 




03546 


34 


05486 


00701 


04040 


D026 


TFM 


L.I.. 


LOOP FOR ADJUSTMENT OF ENTIRE 


03558 


16 


05686 


-0001 


04050 




c 


L»K»» 


MATRIX 


03570 


24 


05686 


11718 


04060 




BNE 


P20 




03582 


47 


03626 


01200 


040 70 


P17 


A 


MATl .NN 




03594 


21 


05491 


05494 


04080 




SK 


DCF 




03606 


34 


05486 


00701 


04090 




B7 


P26-24 




03618 


49 


04094 


00000 


04100 


P20 


TFM 


C0RE,B-9 




03626 


16 


05499 


-7690 


041 10 




RDGN 


DCF 




03638 


36 


05486 


007O0 


04130 




TFL 


Z.ZERO 




03650 


06 


11738 


11713 


041 40 




CM 


K.l 




03662 


14 


11718 


-0001 


U41 50 




BZ 


P23 




03674 


46 


03842 


01200 


04160 


P21 


TFM 


MM.l . . 


COMPUTE Z 


03686 


16 


05676 


-0001 


041 70 




MM 


MM. 10 




03698 


13 
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SAMPLE PROBLEM INPUT 
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^.799382^E-0 ^ 
1 .6565215E-03 

o.oooooooE-qg 

1.6396982E-03 

-1.9032720E-03 
O,O0OO00OE-99 

-1 .6396969E-03 
O,O00O0OOE-99 
O.OOOOOOOE-99 
1.76104^6E-03 
0,000C000E-99 
1.6565215E-03 

-1.7610472E-03 
O.OOOOOOOE-99 

-1.7087824E-03 
0,0000000E-99 
O.OOOOOOOE-09 
1.9032670E-03 
0,0000000E-99 
l,6396g56E-03 

-1.6565215E-03 
^.799382'tE-O^ 

-l,7610433E-03 
O,00OOOO0E-99 
0,OO0OOOOE-99 
1.6396969E-03 
O.OOOOOOOE-99 
1 .7087850E-03 

-1 .761048AE-03 
O.OOOOOOOE-99 

-1.6565215E-03 
O.UOOOOOOE-99 
O.OOOOOOOE-99 
1.7087824E-03^ 
0,O0OOOO0E-99 
1.7610458E-03 

-1.7087850E-03 
0,UOOU00OE-'59 

-1.7610446E-03 
0.OOOOOOOE-99_ 
0,0OO0l5WE-99 



1 .9032720E-03 
O.OOOOOOOE-99 
1 .6396969E-03 
O,0OO00O0E-99 
^.799382'VE-04 

-1 .7610A84E-03 
O.OOOOOOOE-99 

-1.6042594E-03 
1 .7610A72E-03 
O.OOOOOOOE-99 
1.7087624E-03 
0,OOOO00OE-99 
O.OOOOOOOE-99 

-1.70878^0E-53 
0,0000000E-99 

-1.6396956E-03 
1.6565215E-03 
O.OOOOOOOE-99 
1.7610433E-03 
U.UOOOOOUE-99 
0,0000000E-99 

-1.6396969E-03 
O.OOOOOOOE-99 

-1.6565189E-03 
1.7610484E-03 
A. /993824E-04 
1.6565215E-03 
O.OOOOOOOE-99 
0,OOOOOOOE-99 

-1.9032720E-03 
0,OOOOOOOE-99 

-1.6396969E-03 
1.7087850E-03 
O.qoOOOOOE-99 
1.7610446E-03 
O.OOOOOOOE-99 
O,0OOO00OE-99 

-1.7610472E-03 
O.OOOOOOOE-99 

-J.. 70878246-03 

'1.6396969E-03 



O.OOOOOOOE-99 
1.7610484E-03 
O.OOOOOOOE-99 
1 .6042594E-03 

-1.7610472E-03 
O.OOOOOOOE-99 

-1.6565241E-03 
0.C)000000E-99 
4.7993824E-04 
1.7087850E-03 
O.OOOOOOOE-99 
1.6396966E-03 

-1.9032695E-03 
O.OOOOOOOE-'^9 

-1.6396982E-03 
O.OOOOOOOE-99 
O.OOOOOOOE-99 
1.6396969E-03 
O.OOOOOOOE-99 
1.6665189E-03 

-1 .7610446E-03 
O.OOOOOOOE-99 

-1.&04259AE-03 
O.OOOOOOOE-99 
O.OOOOOOOE-99 
1.903:?720E-03 
O.OOOOOOOE-99 
1.6396969E-03 

-1.7C^87g50E-D3 
4.7993824E-04 

-1.7610484E-03 
O.OOOOOOOE-99 
O.OOOOOOOE-99 
1.7610472E-03 
O.OOOOOOOE-99 
1.7087824E-03 

-1.7610458E-03 
O.OOOOOOOE-99 

-1.7087850E-03 
.0000000 E-99 
O.^DTTOOOOE^^ 



1.761047;jh-03 
O.OOOOOOOfc-99 
1.6565241E-03 



O.OOOOOOOE-99 
-1 .7087a50E-03 
O.OOOOOOOE-99 



O.OOOOOOOE-99 
1.7087850E-03 
O.OOOOOOOE-99 



-1 .6565215E-03 
O.OOOOOOOE-99 
-1 .6396982E-03 



1.9032695E-03 
O.OOOOOOOE-99 
1.6396982E-03 

4.7993824E-04 
-1.7610458E-03 



O.OOOOOOOE-99 
1.7610468E-03 
O.OOOOOOOE-99 

-1 .7610446E-03 
O.OOOOOOOE-99 



O.OOOOOOOE-99 - 1 . 65652 1 5E-0 3 



1.7610446E-03 
O.OOOOOOOE-99 
1.6042594E-03 



4.7993824E-04 
1.7087824E-03 
O.OOOOOOOE-99 



O.OOOOOOOE-99 
-1.7087824E-03 
O.OOOOOOOE-99 

1.7087850E-03 
O.OOOOOUOt-99 



-1.9032670E-03 
O.OOOOOOOE-99 
-1.6396956E-03 

O.OOOOOOOE-99 
1.7610472E-03 



1.7610484E-03 O.OOOOOOOE-99 



O.OOOOOOOE-99 
-1.7610472E-03 
O.OOOOOOOE-99 



■1 .6396969E-03 
O.OOOOOOOE-99 
■1.7087850E-03 



1.7610458E-03 
4.799 3824 E-04 
"i. '708785 OE-03 

O.OOOOOOOE-99 
-1.9032695E-03 
O.OOOOOOOE-99 



O.OQOOOOOE-99 
1.9032695E-03 
0.6000bOOE-99 

■1.7087824E-03 
4.799iS24E-04 
-1.7610458E-03 



1.7087824E-03 O.OOOOOOOE-99 
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1.7610433E-03 


O.OOOOOOOE-99 


1.6565215E-03 


O.OOOOOOOE 
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1.7610446E 
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O.OOOOOOOE-99 
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-03 


O.OOOOOOOE 


-99 


1.7087824E-03 


O.OOOOOOOE-99 


1.7610433E-03 










-1.6396969E-03 


O.UOOOOOOE-99 


-1.7087824E-03 


O.OOOOOOOE 


-99 


-1.7610433E 


-03 


O.OOOOOOOE-99 


-1.6565215E-03 


O.OOOOOOOE-99 


-1.7610446E 


-03 


O.OOOOOOOE 


-99 


-1.9032670E-03 


4.7993824E-04 


-1.6396969E-03 


O.OOOOOOOE 


-99 


-1.7087824E 


-03 


O.OOOOOOOE-99 


-1.7610433E-03 


O.OOOOOOOE-99 










0,0000000E-99 


1.6565241E-03 


O.OOOOOOOE-99 


1.6396982E 


-03 


O.OOOOOOOE 


-99 


1.6042594E-03 


0,O00CO00E-99 


1.7610484E-03 


O.OOOOOOOE 


-99 


1.7087850E 


-03 


O.OOOOOOOE-99 


1.6396969E-03 


4.7993824E-04 


1.9032720E 


-03 


O.OOOOOOOE 


-99 


1.7610472E-03 


0,OO0OOOOE-99 


1.6565215E-03 










-1.6565241E-03 


O.OOOOOOOE-99 


-1.6396982E-03 


O.OOOOOOOE 


-99 


-1.6042594E 


-03 


O.OOOOOOOE-99 


-l,7610484E-03 


O.OOOOOOOE-99 


-1.7087850E 


-0 3 


O.OOOOOOOE 


-99 


-1.6396969E-03 


0,00OO0OOE-99 


-1.9032720E-63 


4.7993824E 


-04 


-1.7610472E 


-03 


0,OOOOOOOE-99 


-1.6565215E-03 


O.OOOOOOOE-99 










0,0000000E-99 


1.6396982E-03 


O.OOOOOOOE-99 


1.6565215E 


-03 


O.OOOOOOOE 


-99 


1.6396956E-03 


O.OOOOOOOE-99 


1.7087850E-03 


O.OOOOOOOE 


-99 


1.7610458E 


-03 


O.OOOOOOOE-99 


1.7087824E-03 


O.OOOOOOOE-99 


1.7610472E 


-03 


4.7993824E 


-04 


1.9032695E-03 


O.OOOOOOOE-99 


1.7610446E-03 










-1.6396982E-03 


O.OOOOOOdE-99 


-1.6565215E-03 


O.OOOOOOOE 


-99 


-1.6396956E 


-03 


O,0OO0000E-99 


-1.7087850E-03 


O.OOOOOOOE-99 


-1.7610458E 


-03 


O.OOOOOOOE 


-99 


-1.7087824E-03 


O.OOOOOOOE-99 


-1.7T610472E-03 


O.OOOOOOOE 


-99 


-1.9032695E 


-03 


4.7993824E-04 


-l,7610446E-03 


O.OOOOOOOE-99 










O.OOOOOOOE-99 


1.6042594E-03 


O.OOOOOOOE-99 


1.6396956E 


-03 


O.OOOOOOOE 


-99 


1.6565189E-03 


O.OOOOOOOE-99 


1.6396969E-03 


O.OOOOOOOE 


-99 


1.7087824E 


-03 


O.OOOOOOOE-99 


1.7610433E-03 


O.OOOOOOOE-99 


l.(S.565215"E 


-03 


'O.OOOOOOOE 


-99 


1.7610446E-03 


4.7993824E-04 


1.9032670E-03 










-1 .6042594E-03 


O.OOOOOOOE-99 


-1.6396956E-03 


O.OOOOOOOE 


-99 


-1.6565189E 


-03 


O.OOOOOOOE-49 


-1.6396969E-03 


O.OOOO000E-9?> 


-1.70fl7824E 


-03 


O.OOOOOOOE 


-99 


-1 .7610433E-03 


O.OOOOOOOE-99 


-1.6565215E-03 


O.OOOOOOOE 


-99 


-1.7610446E 


-03 


0,O0O000OE-99 


-1.90 32670E-03 


4.7993824E-04 










l.OOOOOOOE 00 


O.OOOOOOOE-99 


l.OOOOOOOE 00 


O.OOOOOOOE 


-99 


l.OOOOOOOE 


00 


O.OOOOOOOE-99 


l.OOOOOOOE 00 


O.OOOOOOOE-99 


l.OOOOOOOE 


00 


O.OOOOOOOE 


-99 


l.OOOOOOOE 00 


a.OOOOOOOE-99 


l.OOOOOOOE 00 


O.OOOOOOOE 


-99 


l.OOOOOOOE 


00 




l.OOOOOOOE 00 


O.OOOOOOOE-99 


l.OOOOOOOE 


00 


O.OOOOOOOE 


-99 
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SAMPLE PROBLEM OUTPUT 
RO^ 01 



1»^66 3088E+03 


6.45 75 44 7E + 02 


-3^7460682E+02 


7*4^sn'if»PF + n 1 


_ 1 




-1 •C332652E + 02 


— 3»7460671E + 0? 


7 • 46 49988E + 01 


1 • / 0\J^£,L.~\j£. 


_ 2 




-6«6636106F+01 


-1 •4405735F+02 


— 1 .741 S'^P^F + np 




~~6 






-?.4 905A??F+ni 


— T.?7'iA1P P F ■♦■0 P 












KUW 0<i 








-6«4575471E+02 


1 •466308 1 E + 3 


— 7»4650 132E + 01 


— . 7^ AH A 7 AP + 9 




n 5 A n F * n 9 




— 7.4A497PlF+ni 


— ^» 74AAAAAF*AP 
j% 1 HO u DO DC TU A 






n 7 B 1 1 P A A 5 


1 • A 4 5 7 1 7 F + 2 


— 6»6636169E+01 


1 •OTSPAnPF+flP 


— 1 - 7A1«^'^1 7F4.n9 






— 6»663607AE + '^1 


1 .P7fiA1A ^ F + n 5 


— P-4QD^41 OF*ni 












ROW 03 












1 •a < i IWUt+UJ 


!)»:>15oOoofc + 02 


— 3 


• 7460658E+02 


7 ^ A n "^ri Qp 4.n 1 


— 1 Q ft 4 ? ^ ^ w r J. fi •> 


— 1 7 ft P ZlF A n 1 




9 


• Vooo iOOt — 1 






«.A-AA'^A4n4F*rtl 


— 1 _ A AO 1^ 7 "3 F * n "5 




C« A ^ *T C C .4. A T 


~1 «4754309F40? 


— A A A A 4 n P J. 1 


— 1 440 ^ 7 P AF J. r\ P 












ROW 04 








•"7»4650461E + 01 


— 3»7460686F + 02 


— 'i.'^l'ift1'^7F + np 






. A A A P^ A 1 




Q i7ft'aiAACj.ni 


— 1 Oft/i7CROF^rtO 

""I •704i'PPC3t + U^ 




-2 


•4531696E+02 


9» 178 3322E + 01 


-1 •9842512E+02 


1 .44n5740F+n? 






-.Zi7*iA'5A7FAA0 


-5»2650584E+01 




— A.AA'^SftA4F4-01 












ROW 05 








- 1 • 74 1 5 2 92E+02 


-1 •03 52 6 1 9E+02 


— '?.74A0AA9F + 0P 








6 • 45 75 5 1 8E + 02 


— fe.ftft'^Ani iF+ni 


— 1 .440^ 7A4F + np 


— 1 - Q 7 70Qft F 7 


_ 2 




— J» /4oUo4ot+U^ 


f •4opU3 /ot+Ul 


— 2»4905406E + 01 


-1 •2 756 13^E+02 


-6 


•6636140E+01 


-1 •4405719E + 02 


-1 •7415260E + 02 


— 1 i03'iP'591F4-0P 












RDW Oft 








1 •0352616E + 02 


-1 •7415296E+02 


-7^4650439E+01 




"6 




1 - 4 AA ^ Ofl ft F -4.n 




— A^AA'SAn'^OFj.n 1 
O.OOJDU57C+U1 




"~ 1 




1 — 7 « 465 1 5 4E + 1 


— '^-74AnA^nF4-np 


1 - P7<^A1 "^ftF^np 

A»A l-?OAPOU^^A 




1 




A A ii Q A n "7 rt c ^ n 1 


1 n^ROK71FAr\"l 


-!• (4l5^ppt + U<i 












ROW 7 








-3.7460671E+02 


7«4649902E+0l 


— 1 * Qfl 4 P S AF 4.np 


— 0-1 7R "i^ 5 Q^F 4.h 1 


-* A 




-1 •4405726E + 02 


1.571 1990E+03 


5,3258069E+02 




9 




-5.2650722E+01 


— 1 •^754312E+02 


— 3 • 7460684 E + 02 








-9» 1783250E + 01 


-6»6636270E+0l 


-1^4405735E+02 




















-7.4650069E+01 


-3.7A60695E+02 


9.1783092E+01 


-U9842557E+02 


1 


.4405716E+02 


-6.6636439E+01 


-5,3158119E+02 


1.3711986E+03 


-9,99«.7&l0e-01 


-2 


.4531720E+0? 


1 .4754326E + 02 


-5,2650556E+01 


-7.^649847E+01 


-3.7460674E+02 


9 


.1783371E+01 


-1 ,9842532E + 02 


1 «4405762E+02 


-6.6635987E+01 












ROW 09 




_____ 




-1 .~90780^f3E+02 


-i .0279786E+02 


-2.^531705E+02 


'9.9865530E-dl 




.9078013E+02 


-1.0279773E+02 


-2.4531715E+02 


9»9871160E-01 


1.7090600E+03 




.5992369E+02 


-2.4531696F+02 


9t9989500E-0] 


"lt907Q038E+Q^ 


-liO?79791E+0? 


-?. 


.453169QE+02 


9.9856500E-01 


-1 •9078006E+02 


-1.0279789E+02 
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1 •0279799E + 02 
•1 .9Q78006E-t-02 



■1 .9078048E+02 



ROW 10 
-9^9853900E-01 
-2t45317Q9Ef0? 



-2.4531700E+02 
- 4t 5992 4 12 Et0 2 



1.0279769E+02 
l« 7Q 9Q599EtQ3 



-9^9883000E-0l 
-2^4531695E+02 

-6^6636208E+01 
7^4650295E+01 
1 .571 1988E-f03 



-2.4531697E+02 
1 ,0279788E+02 

-1.4405735E+02 
-5,2650751E+01 
^t31 ^81 ^2E-^02 



1^0279773E+02 
-1^9078012E+02 
ROJfl/ 1L__„ . 
-1^9842539E+02 
-1.4754308E+02 
- »f^ 6 35904E -»- 1 



•1 •9078034E+02 -9 • 9844900E-01 



-9^1783212E+01 
-2.4531727E+02 
•1 .4A05770E»02 



-3.7460651E+02 
9.9929000E-01 
•lt98425Q2E-*-Q2 



-9.1783235E+01 -3 .7460622E+02 



1 .4405748E+02 
-3^7460648E+02 
-5^3158175E+02 



-6.6636059F+01 
l,4754315E+02 
1.5711989E+03 



7^4650304E+01 

ROW 12 
9,1782980E+01 
-5^2650645E+01 
1.4405751E+02 



-1,9842528E+Q2 
-9,9896500E-01 
-6^6636000E+01 



-7,4650240E+01 
-2.4531733E+02 
9.1783167E+01 







ROW 13 












-1.7415305E+02 


-1 .0352627E+02 


-6^6636157E+01 


-1 


.4405742E+02 


-2 


.4905349E+01 




-1^2756150E+02 


-3.7460688E+02 


7^4650094E+01 


-1 


.9078033E+02 


-1 


.0279804E+02 




-6^6636266E+01 


-U4405721E+02 


1 .4663086E + 03 


6 


.4575474E+02 


-3 


.7460701E+02 




7,4650331E+01 


-1 •7415290E+02 


-U0352640E+02 
ROW 14 













1 •0352633E + 02 
-2^4905331E+01 

1 •4405739E+02 
-3,7460682E+02 



-1 •7415331E + 02 
-7.4650053E+01 
-6^6636280E+01 
K0352623E + 02 



-6^6636448E+01 -1 • 4405 744E+02 



1^4405762E+02 
-3.7460692E+02 
-6.4575497E+02 
-1^7415284E+02 

ROW 15 
■5^2650617E-»-01 



-6.6636246E+01 
1.0279782E+02 
1.4663087E+03 



1.2756160E+02 
■1.9078025E+02 
-7.4650269E+01 



-1.4754318E-t-02 -6.6636081E-»Q1 



-1 .4405719E+02 
•1 •9842514F+02 
5.3158078E+02 



-1 .9842575E + 02 
-9, 1783117E + 01 
-3.7460642E+02 



1.4405741E+02 -6 .6635963E+0 1 



-9.1782889E+01 
-3^7460697E+02 

7^4650187E+01 
ROW 16 

1.4754329E+02 



-2.4531698E+02 
7.4650370E+01 



9,9858700E-01 
1.5711991E+03 



•6.2650646E+01 1 .4405728E+02 



9.1783044E+01 


-U9842511E + 02 


-7^4650249E+01 


-3 


•7460705E+02 


-5 


.3158132E+02 


1.57] 1989E+03 


-7,465018CE+01 


-3^7460667E+02 
ROW 17 










-2^4905246E+01 


-1 ,2756135E + 02 


-6^6635970E+01 


-1 


•4405744E+02 


-1 


.7415252E+02 


-1 •0352599E + 02 


-6.6636280F+01 


-1.4405766E+02 


-1 


.9078030E+02 


-1 


.0279773E+02 


-3^7460657E+02 


7^4650430E+01 


-1.7415299E+02 


-1 


•0352592E+02 


-3 


.7460661E+02 


7^4650420E+01 


1 •4663090E + 03 


6^4575500E+02 
ROW 18 










1 .2756157E + 02 


-2^4905077E+01 


1 •4405743E + 02 


-6 


•6635999E+01 


1 


.0352625E+02 


-1 • 741 5264E + 02 


1 •4405735E + 02 


-6^6636260E+01 


1 


•0279760E+02 


-1 


.9078034E+02 


-7^4650438E+01 


-3^7460659E+02 


1 .0352623E + 02 


-1 


•7415307E+02 


-7 


.4650306E+01 


-3.7460659E+02 


-6^4575543E+02 


1 .4663090E+03 











2 5 



VARIABLE 



SOLUT ION 



01 1.9830998E+01 

02 -6.95283^0E+01 

03 -6.1050^00E+00 

04 -6.2656100E+01 

05 1.9832610E+01 
Of! -6.9529790E-t-01 

07 -6.1058000E+00 

08 -6.2655270E+01 

09 -3.5329000E+01 

10 -5.2727110E+01 

11 -6. 10«»4400E+00 

12 -6.2656883E-t-01 

13 1 .9830690E + 01 

lA -6.9528290E+01 

15 -6.10522O0E+O0 

16 -6.265605DE+01 

17 1.9832600E+01 
18 -6»9529730E-»-01 
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